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We use a continued fraction expansion of the sign-function in order to obtain 
a five dimensional formulation of the overlap lattice Dirac operator. Within 
this formulation the inverse of the overlap operator can be calculated by a 
single Krylov space method and nested conjugate gradient procedures are 
avoided. We point out that the five dimensional linear system can be made 
well conditioned using equivalence transformations on the continued fractions. 

1 Introduction 

Let us start with noting the overlap Dirac operator D describing chirally 
symmetric fermions on the lattice pQ, 



where H(—m) is a hermitian Dirac operator with a negative mass parameter 
— m of the order of the cut-off of the lattice theory, m ~ 0(a). A bare quark 
mass [i is most conveniently introduced as D(/i) = (1 — \i)D + In order to 
calculate efficiently the sign-function in eq. Q one can use a rational approx- 
imation sign(a:) ~ R n ^ m (x) where R n , m is a (nondegenerate and irreducible) 
rational function with algebraic polynomials of order n and m as numera- 
tor and denominator, respectively. Rational approximations usually converge 
much faster with their degree than polynomial approximations, but of course 
for our specific application it might still be much more expensive to apply the 
low degree denominator of the rational function than to apply a high order 
polynomial. However, noting that a rational function can be written as a par- 
tial fraction by matching poles and residues, i.e., R n , m {x) ~ x Ck/(x 2 +qk), 
one can use a multi-shift linear system solver where the convergence is gov- 
erned by the smallest of the shifts qu- The overall cost is therefore roughly 
equivalent to one inversion of H 2 . 
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In order to do physics we need to compute inverses of the overlap operator 
D(ijl) to obtain propagators, fermionic forces for Hybrid Monte Carlo, etc. If 
we consider the multi-shift linear system above, we realise that the inversion 
of D(jjl) leads to a two-level nested linear system solution, which is rather 
cumbersome and forbidding. It is well known how this can be avoided [30]: by 
introducing a continued fraction representation of the rational approximation 
and auxiliary fields, the non-linear system D(p)ip = \ can be unfolded into 
a set of systems linear in H . The auxiliary fields can be interpreted as fields 
living in a fictitious fifth dimension and in this way the nested 4D Krylov 
space problem reduces to finding a solution in a single 5D Krylov space. One 
can also regard the auxiliary fields as additional fermion flavours which, when 
integrated out, generate an effective Dirac operator equivalent to -D(/i). 



2 Matrix representation of the sign-function 

Consider a rational approximation to the sign-function and expand it as a 
continued fraction 3 , 
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If we rewrite the linear system 
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using appropriate auxiliary fields 
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i>2 n we obtain the system 
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By performing explicitly a Schur decomposition of the matrix it is easy to see 
that the inverse of eq.J5J is just the (1, l)-component of the Schur complement. 

The rational function can also be mapped into a so-called simple continued 
fraction 



Note that since polynomial approximations can also be expanded into continued 
fractions all our considerations apply to them as well. 
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R2n+l,2n{x) = fax 



which leads to a different matrix 
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and we find that the ctfj's and fa's are related through 

fa = "0! fa = ) ■■ ■ ,fa — 3 7 ■ 
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In order to understand the relation between different continued fraction rep- 
resentations of the same rational function in detail we need to take a closer 
look at the properties of continued fractions (see e.g.^JE])- 



3 Continued fractions 

A generic (truncated) continued fraction ^ is conveniently written as 
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Simple continued fractions have the property ai = l,i 
usually writes 

-7T- = [fa] faifa,---: fa] ■ 



1,. 



and one 
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Continued fractions are widely used in many areas of physics and mathe- 
matics, in particular also in number theory Finite continued fractions provide 
an alternative representation of rational numbers and form the basis of ra- 
tional approximation theory. Infinite continued fractions on the other hand 
can be used to represent irrational numbers. Some numbers have beautiful 
continued fraction expansions while others have very mysterious ones. Let us 
quickly note a few examples for our amusement: 
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0= [1; 1,1, 1,1,...], 
V2= [1;2,2,2,2,...], 
e= [2;1, 2,1, 1,4, 1,1,6,1, 1,8,1, 1,10,1,...], 
tt= [3; 7, 15,1, 292, 1,1,1,2,1,3,1,...], 

where (/> = |(1 + \/5) is the golden mean. It is interesting to note that there 
is no regular pattern known for 7r, and it is not known why this is so. 

Evaluation of continued fractions can be done through forward or back- 
ward recurrence algorithms, and the former makes use of the intimate relation 
between continued fractions and the coupled two term relations 

A-x = 1, A = fa, £_! = 0, B = 1, (11) 
Ai = fcAi-! + anAi-2, i = 1,2,3,..., (12) 
Bi = fiiBi-i + aiBi-2, i = 1,2,3,..., (13) 

which can equivalently be written as an iterative matrix equation, 

AA (A { _ x A^ 2 \ ( pi 
Bi J I Bj_i Bi_ 2 J { ai 



(14) 



There is also an interesting connection between continued fractions, Moebius 
transformations and the corresponding unimodular matrices. 

The natural arithmetic operation for continued fractions is inversion and 
the corresponding rule is particularly simple: 

{R . R r i _ / [0;/3 ,/3i,...] if A) 7^0, . . 

l#>.0i>-"] -{[/3 i; /3 2 ,...] if/3 = 0. (15) 

It is also helpful to write down the rule for the multiplication of a continued 
fraction by a constant c, 

c- [/3 ;/3i,...,/3„] = [c-/3 ;-,c-/3 2 ,^,c-/3 4 ,...]. (16) 

c c 

Most important for our purpose, however, is the equivalence transformation 
of a continued fraction which is stated in the following theorem: 

Theorem 1. Two continued fractions Po + j^ + ■■ - + an d f3' o + j0r + . . . + ^r 

are equivalent iff there exists a sequence of non-zero constants c n with Cq = 1 
such that 

a'i = CiCi-icti, i = 1,2,3, . . . ,n, (17) 
/?•= Cift, t = 0,l,2,...,n. (18) 

The theorem is easily seen to hold true by explicitly writing out the full 
continued fraction, 
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Po + = 0o + — (19) 

OL2 - CiC 2 «2 
Pi + *01 + 

Pi + C 2 P 2 + 

^ a n c n -\c n a n 

Pri CnPn 

and it is also clear that in terms of approximants we simply have ^ = 

C1C2 c"a" ' ^° we ^ n< ^ a §i ven rational function corresponds to an equiv- 
alence class of continued fractions and the class is parametrised by the (non- 
zero) coefficients Cj. While the equivalence transformation itself appears to be 
rather trivial, and indeed leaves the value of the continued fraction invariant, 
it affects the spectrum of the corresponding matrix representation [3] . 

In the analytic theory of continued fractions one considers continued frac- 
tions of the form 

ftM + ^ + ^ + "-. ( 2 °) 

\Pl{z) \P2{z) 

i.e., the coefficients are functions of a complex variable z. So-called J-fractions 
are of the special form 

[nz + si,r 2 z + s 2 ,r 3 z + s 3 , . . .], (21) 

where Ti 7 Si are complex numbers with 7^ 0. One can show that the n-th 
approximant of a J-fraction is an element of R n -i^ n . Conversely, let 

Pn-i(z) = Pl z n - 1 +p 2 z n - 2 + ...+p n , (22) 

Q n (z) = q z n + q 1 z n - 1 + ...+q n , (23) 

then we have 

-TTT^T = I r i z + s ii r i z + 8 2 ,.. .,r n z + s n ] (24) 

Qn\Z) 

where r^, Si are uniquely determined by Pi,qi- We can now replace z by —z, 
apply an appropriate equivalence transformation and, using the fact that 
P n _x(z)/Q n (z) is odd, we find Sj = from uniqueness. Therefore we can 
always write 

P„_i(z 2 ) 

z 7\ i ?\ =[hz,k 2 z,k 3 z,...,k m z], (25) 

Qn(Z 2 ) 

where m = In — 1 or m = 2n according as Q n (0) = or Q„(0) =^ 0. 



4 Application to the overlap operator 

We are now in a position to apply our knowledge to find the solution to 
the equation j— 75_D(/i)^J = X- Collecting the results from the previous two 
sections we obtain an equivalence class of five dimensional linear systems 
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where the fcj's and n are uniquely determined by the given rational approxima- 
tion to the sign-function, the q's parametrise the corresponding equivalence 
class and A = j^j;- It is now crucial to see how the spectrum of the five di- 
mensional matrix depends on the free parameters Cj as we already emphasised 
in the last section. While for a generic set of parameters the five dimensional 
system is usually ill-defined, the condition number can be kept under control 
with a clever choice of Cj's jSlEj enabling one to optimise the matrix, e.g., for 
fast inversions. As was pointed out in |3] the equivalence transformations can 
be understood as a block Jacobi preconditioning without any computational 
overhead. The particularly simple structure of the five dimensional operator 
allows improvements in various directions: one can easily change and optimise 
the hcrmitian overlap kernel H or apply various well known preconditioning 
techniques such as an even-odd or ILU decomposition [0]. 

It is instructive to see that the first auxiliary field disentangles 75 from the 
sign-function, i.e., (A ■ 75 + sign(H))ip — \ maps into 



A- lb 



-c 2 sign(i/) 



(26) 



where we used sign = sign - . Additional fields are then only used to gener- 
ate the sign-function (or its approximation, respectively). We now have two 
systems which can be easily solved, 



yielding recursion relations for ip and (j>, 
V (i+1) = ^75(x-sign(ff)V> (i) ), <M i+1) 



1 



sign(H)ip, 



^■sign(iJ)7 5 (ix- 



(27) 



(28) 



Equivalently, there are recursion relations for the residuals and one can show 
that 



'■0,01 



V1 + /J 
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Of course one can use a similar recursive scheme for the case where sign(_ff ) 
is expressed as a continued fraction and one has several auxiliary fields. 

Projection of eigenvectors of H close to is a valuable tool to improve 
approximations to the overlap operator and here it is straightforward to im- 
plement. However, we wish to point out that in this formulation it might not 
be necessary at all. Consider the linear system in the lower right corner, 
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C2n~lC2n4>2n-l ~ cl n k 2n H (f>2 n = 0. (30) 

For a typical rational approximation to the sign-function we have k 2n 3> 1 
and with the choice C2 n — &2n+i "C 1. So it turns out that the system 
in eq.ljSOJ is essentially equivalent to finding eigenvectors of H close to zero, 
i.e., the two Krylov spaces possibly have a large overlap. Indeed a trunca- 
tion in the fifth dimension, i.e., of a given continued fraction, changes the 
approximation only in the neighbourhood around H = and therefore pro- 
vides a natural truncation scheme for approximations of fixed accuracy. It also 
opens up the possibility of applying successively better approximations to the 
overlap operator which are ultra-local in five dimensions. 



5 Summary and outlook 

We have shown how to use a continued fraction expansion of the sign-function 
in order to obtain a five dimensional formulation of the overlap lattice Dirac 
operator. We have pointed out that the operator can be made well condi- 
tioned using equivalence transformations on the continued fractions. It is now 
important to investigate in detail strategies to exploit this freedom in prac- 
tical applications and such a study is under way [S]. If successful, and first 
results indicate that this is indeed the case, the formulation would provide a 
valuable alternative for the simulation of dynamical chiral fermions. 
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